Roton softening and supersolidity in Rb spinor condensates 
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Superfluids with a tendency towards periodic crystalline order have both a phonon and roton like 
spectrum of collective modes. The softening of the roton spectrum provides one route to a supersolid. 
We show that roton softening occurs in 87 Rb spinor condensates once dipolar interactions and spin 
dynamics are taken into account. By including the effects of a quasi-two-dimensional geometry and 
rapid Larmor precession, we show a dynamical instability develops in the collective mode spectrum 
at finite wavevectors. We construct phase diagrams showing a variety of instabilities as a function of 
the direction of the magnetic field and strength of the quadratic Zeeman shift. Our results provide 
a possible explanation of current experiments in the Berkeley group Phys. Rev. Lett. 100:170403 
(2008). 



The experimentally elusive supersolid state has long 
been of great interest as an exotic quantum state of 
matter. Such a state offers the possibility of periodic 
crystalline order breaking translation invariance famil- 
iar from solids and long range phase coherence breaking 
global gauge invariance familiar from superfluids coex- 
isting in the same material [j], 0, H, 0, [H, El • Theoretical 
interest in the supersolid state has a long history dat- 
ing back to early studies of superfluidity in 4 He. Landau 
suggested that the excitation spectrum consists of two 
parts: a sound like long wavelength mode and a roton 
spectrum e(p) = A + (p — po) 2 /2m for p ss po indicative 
of a tendency towards crystalline order. In a system with 
both a phonon and roton spectrum, the softening of the 
roton gap where A approaches zero provides one means 
of possibly realizing a supersolid 0,11] ■ Recent experi- 
ments suggesting the existence of the supersolid phase in 
4 He are a subject of intense debate 0, [l(| [TJ, • 

In this paper, we demonstrate roton softening takes 
place in quasi-two-dimensional F = 1 ferromangetic con- 
densates such as ultracold 87 Rb once dipolar interactions 
and spin dynamics are taken into account. Pattern for- 
mation due to dipolar interactions is well known in con- 
densed matter physics [13|, but spinor condensates such 
as 87 Rb present several novel effects without condensed 
matter analogs. In addition to dipolar interactions, the 
spin degrees of freedom also experience competing in- 
teractions including spin dependent contact interactions 
and the quadratic Zeeman shift [14j. Moreover, the spin 
dynamics of rapid Larmor precession and confinement 
to a two-dimensional geometry play an important role in 
modifying the effective dipolar interaction. Several works 
have previously analyzed the role o f dip olar interactions 



in polar molecules [15|, [16|, ll7|, [18|, [19( and spinor con 



densates [2(J [2l[ [22[ . However, these studies do not take 
into account the combined effects of dynamical spin de- 
grees of freedom, rapid Larmor precession, and reduced 
dimensionality. 

We show that for an initially uniform ferromagnet, the 
excitation spectrum has roton like parts for both spin and 
density branches (generally each branch invloves both 



spin and density degress of freedom, so we define them 
by their behavior in the long wavelength limit) with the 
roton on the spin branch becoming imaginary at finite 
wavevectors. Macroscopic occupation of such roton exci- 
tations should lead to a state which breaks global gauge 
invariance, spin rotational and translational symmetries. 
Hence such instability suggests the likely formation of 
a supersolid phase. This instability has a simple phys- 
ical origin in terms of lowering the classical dipolar in- 
teraction energy of a uniform ferromagnet through pe- 
riodic modulation of the magnetization F as shown in 
Fig. Q] However, the modulation direction is differ- 
ent for the longitudinal component along the direction 
of the magnetic field B versus the perpendicular compo- 
nents because only the latter precesses around the ex- 
ternal magnetic field (see discussion below). Different 
components of the spin are conjugate variables which 
couple to each other through the commutation relation 
[F 1 ,^] = ie v F , a rich variety of instabilities leading 
to possible striped and checkerboard supersolid phases 
can arise. This simple physical picture also gives a rough 
estimate for the length scale A of the instability by equat- 
ing the kinetic energy cost with the gain in dipolar in- 
teraction energy l/2mA 2 + \/2md 2 = 2irgdno/3 where m 
is the mass, gd gives the strength of dipolar interactions, 
no is the density, and d n is the thickness of the conden- 
sate. With typical experimental parameters d — 2pm 
and gdno = 10 Hz this gives an estimate of A sa 10/i 
m. The length scale and structure of unstable modes 
agree quantitatively with current experiments on dipolar 
effects in spinor condensates [231 ] . 



RESULTS 

We consider a quasi-two-dimensional spinor conden- 
sate as shown in Fig. [2] The unit vectors x, y are in the 
plane while n is out of the plane. A uniform magnetic 
field points along B in the n, x plane at an angle a with 
respect to n. The initial condensate is either prepared 
with a uniform magnetization or with a non-uniform spi- 
ral spin texture with magnetization winding along k and 
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FIG. 1: We consider a quasi- two-dimensional spinor conden- 
sate with magnetic field B along the 3 direction in the plane, 
and uniform magnetization F. Larmor precession of F about 
B means F spends half the time along the 1 direction in the 
plane and half the time along the 2 direction out of the plane. 
Magnetization fluctuations 8F with periodic modulation can 
lower the effective dipolar interaction energy and drive an 
instability towards a possible supersolid state. 5F along 3 
(parallel to B) favor modulation along 1 (perpendicular to 
B). In contrast, 5F along 1 or 2 (perpendicular to B) favor 
modulation along 3 (parallel to B). All components couple to 
each other through canonical commutation relations and give 
rise to a variety of instabilities towards striped phases modu- 
lated along 1 or 3 as well as checkboard phases modudulated 
along I and 3. 



wavevector |«|. The Hamiltonian is given by 



H = I d 3 xVt 

d 3 x 



-—-fi + B B-F + q[B-F 
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+ J d 3 xd 3 x'^ D (x~x') : ^j 5 ****^'** : 

(1) 

where : : denotes normal ordering. We denote with 
a = 1,2,3 as annhilation operators for F = 1 bosons 
with mass m and F hyperfine spin operators with F l - k = 
—ieijk- Throughout, we use a matrix notation with sup- 
pressed indices where *, T, and f denote the complex 
conjugate, transpose, and the conjugate transpose, re- 
spectively. For example, \& is a column (row) vector 
while F % is a matrix. 

The chemical potential /i is a Lagrange multiplier con- 
trolling the density n^rj = (M/tlt) and we work with fixed 
longitudinal magnetization tlzdIb = (9f'B ■ F*&). The 
magnetic field induces Larmor precession about B at a 
frequency Bq and a quadratic Zeeman shift q. With typ- 
ical magnetic fields B of zero up to hundreds of mG, 
q = 70 Hz G -2 B 2 LJ] ranges from zero to tens of Hz. 
AC stark shifts can further tune q, in particular to nega- 
tive values. A harmonic trapping potential along n con- 
fines the condensate to a thickness d n . We take typical 
values of Bq/2it = 115 kHz and d n — 2 /im [23]. 

The spin independent and spin dependent contact in- 
teraction strengths are given by go — 4irh 2 ao/m, g s = 
Airh 2 (ao — 02) /3m [24] in terms of the s-wave scatter- 
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FIG. 2: Schematic of experiment with unit vectors x, y (n) 
in the plane (out of the plane). The magnetic field B (spiral 
wavevector k makes an angle a (/3) with respect to n (x). 
The effective dipolar interaction takes into account Larmor 
precession of the magnetization F and the confinement of the 
condensate along b to thickness d n . 



ing lengths cif for two atoms colliding with total an- 
gular momentum F. For 87 Rb, ao = 101. 8ob and 
<i2 = 100. 4a s where is the Bohr radius [25| giving 
positive g s and ferromagnetic interactions. The dipolar 
interaction strength is given by g^ — ^ag 2 F ^ 2 B where /io is 
the vacuum permeability, gp is the Lande g-factor, and 
Hb is the Bohr magneton. The dipolar interaction tensor 
is given by 

4-7T r 
"3 



h£ D (x 
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3i l i J 



hUk) 



6 lJ 



3k l k J 
(2) 



in real and momentum space with the Fourier transform 
regularized as in Methods: Dipolar Interaction. For typi- 
cal peak three-dimensional densities of n^u = 2.2 x 10 14 
cm -3 the interaction strengths are gpn^f) = 1.7 kHz, 
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g s n 3 D = 8 Hz, and g d n 3D = 10 Hz 

Notice the clear separation of energy scales for the 
above Hamiltonian. The quadratic Zeeman energy, spin 
dependent contact interaction, and dipolar interaction all 
compete at the lowest energies on the order of tens of Hz. 
The spin independent contact interaction and harmonic 
trapping along h are both hundreds of Hz implying the 
density and out of plane dynamics are effectively frozen. 
At the highest energies, Larmor precession occurs at hun- 
dreds of kHz implying one should average over this rapid 
precession for the low energy properties. 



Dipolar Interaction 

As illustrated in Fig. O the bare dipolar interaction 
is modified by spatial and time averaging due to confine- 
ment along n and Larmor precession, respectively. Here 
we briefly summarize the resulting effective dipolar inter- 
action that emerges with more details in Methods:Dipolar 
Interaction. 

We first consider the effect of confinement along h. 
Since the thickness of the condensate d n and the spin 
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healing length £ are comparable d n , £ ~ 2 /im, we assume 
the condensate is frozen along n and take 



(3) 



where on the left-hand (right-hand) side is a three- 
dimensional (two-dimensional) field. Here p{x n ) is the 
normalized to J dx n p(x n ) = f and for dcfmitcncss we 
take a gaussian form p(x n ) = exp(— x^/2d^)/ y / 27rd| 
with cc„ the coordinate along h. Integrating over x n 
in Eq. [T] gives a two-dimensional Hamiltonian with 
d 3 x -> ePx, .go -> g o p(0)C, g s -> g s p(0)C, g d -> gd,p{0)C 
where the constant C = l/v2 is determined by nor- 
malization. However, the dipolar interaction tensor 
h l ^ D (Ax) — > h l 2 D (Ax) is modified appreciably. 

In addition to confinement along n, rapid Larmor pre- 
cession also modifies the dipolar interaction tensor. In 
experiments, the strong uniform component of the mag- 
netic field causes precession of magnetization along B 
while a weak gradient along k induces a spiral order of 
the magnetization [23| . We take this into account by go- 
ing to a co-moving frame via the unitary transformation 



R(t, x)^ x , R(t, x) = exp iO(t, x)B ■ F (4) 



with 6(t,x) = —Bot + \k\k ■ x where Bq is the uni- 
form component of the magnetic field and |«| is the 
spiral wavevector. Applying this unitary transforma- 
tion simply yields the substitutions Bq — > — i\K\/mk ■ V, 
q — > q + \ k\ 2 /2m which arise from the Berry's phase and 
kinetic energy terms. Time averaging over the rapid pre- 
cession then yields a modified dipolar interaction tensor 



We focus on the resulting two-dimensional time- 
averaged interaction in momentum space given by 



(5) 



with the momentum independent projection operators 



Pg = B l &, 



p± = s ij - ir ii' (6) 



and recall FK — —ieijk while the functions hfuik), 

h>2£)(k), h,2 D {k) carry the momentum dependence and are 
given explicitly in Eqs. [2H HS1 ESI of Methods: Dipolar In- 
teraction. Notice the spin dependent part carrying the 
i,j indices only depends on B and not k. Essentially, 
time-averaging over the fast Larmor precession selects a 
preferred direction in spin space along B. Compare this 
to the bare dipolar interaction of Eq. [21 where the spin 
and momentum dependence are not separable. In real 
space, the bare dipolar interaction (without averaging 
over Larmor precession) favors spins aligned head-to-tail 
and anti-aligned side-by-side. In particular, small distor- 
tions of a mean-field condensate with uniform magnetiza- 
tion in the plane are energetically unfavorable since they 
destroy the favorable head-to-tail order already present. 



These considerations change once rapid Larmor preces- 
sion is included. Precession causes the energetically fa- 
vorable head-to-tail order to rotate into the energetically 
unfavorable side-by-side order half of the time. Heuristi- 
cally, this gives rise the instabilities of the effective dipo- 
lar interaction. Fig. [1] illustrates the fluctuations SF 
for the magnetization F that lower the effective dipo- 
lar energy with B in the plane. SF parallel to B favors 
modulation perpendicular to B. This is as expected even 
for the bare dipolar interaction which favors head-to-tail 
alignment (uniform parallel to B) and side-by-side anti- 
alignment (modulation perpendicular to B). However, 
both components of 5F perpendicular to B favor modu- 
lation parallel to B. For SF in the plane and perpendic- 
ular to B, this can be understood in terms of the bare 
dipolar interaction favoring head-to-tail alignment (uni- 
form perpendicular to B) and side-by-side anti-alignment 
(modulation parallel to B). Larmor precession affects SF 
out of the plane and perpendicular to B the most as pre- 
cession rotates this component into the plane half the 
time. This leads to the same type of behavior as SF in 
the plane and perpendicular to B: uniform perpendicu- 
lar to B and modulation parallel to B. Compare this to 
arguments appealing to the bare dipolar interaction for 
this component which suggest modulation along both di- 
rections since SF out of the plane gives the energetically 
unfavorable aligned side- by-side arrangement. 



Collective Modes 

We now turn to quantitative analysis of the collective 
mode spectrum focusing on the two-dimensional time- 
averaged case. Starting from mean field solutions with 
uniform magnetization, we study the collective mode 
spectrum describing its small fluctuations. For details 
see Methods: Collective Modes. Recall we transformed to 
a frame co-moving with possible spiral order so that co- 
moving frame uniform states describe both lab frame uni- 
form and spiral states. 

We take B as the quantization axis and parametrize 



sin(fe) 



■ ie im+i» g cos (^ + lX3 ) 

ie ir, S +iu 3 sm((fe + lX3) sin(pa) 
-y/cosh(2Xx 

e 1 ^ cos(p s ) 



(7) 



with n the two-dimensional density, rj the global phase, 
p, Xi v controlling the magnitude of the magnetization, 
and 4> the orientation of the transverse magnetization. 
We take tyg = \& independent of x and show in Meth- 
ods: Collective Modes these mean- field states only depend 
on 

9± =9s-9dhi D {0) (8) 



2.9±«-3dC" 
where we use n%D — np(0). 
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To study collective modes, we take = * + S&g 
and find the linearized equations of motion for the fluc- 
tuations <5*£ given by 



id. 



' 6 *k ' 




[ M % 








— fe_ 










J** - 

— K. 



(9) 



with Afg and iVg given in Methods: Collective Modes. The 
ansatz 5^f^(t) ~ e l "£* gives an eigenvalue equation for 
the excitation energies wr. Analysis of gives the spec- 
trum for small fluctuations above the mean field solution. 
In particular, imaginary lu^ indicates a dynamical insta- 
bility where such fluctuations grow exponentially. 

The above collective mode analysis describes the gen- 
eral case for arbitrary fs, \k\ ^ 0. Here we focus on the 
analytically tractable case of zero longitudinal magneti- 
zation and no spiral order Jb, \k\ = and later discuss 
the effects of finite /b, \k\. Instead of the three compo- 
nent complex field tyg, it will be convenient to use the 
six component real field 



$2 = [ns rjg p s v 3 <fo xs] ' 



(10) 



where Eq. [7] defines as a function of For the 
collective mode analysis, we take 



<9* 
5* 



(11) 



with the derivative evaluated at the mean-field param- 
eters and derive equations of motion for <53?£ from Eq. 
The (p and x modes decouple from the n, rj, p, and v 
modes. Since <f> controls the orientation of the transverse 
magnetization and x controls the magnitude of the lon- 
gitudinal magnetization, they form a conjugate pair of 
variables which we denote as the spin mode. The other 
four degrees of freedom form two pairs of conjugate vari- 
ables which we denote as the charge and magnetization 
modes. 

At long wavelengths, we find for the spin, charge, and 
magnetization modes 



= J[1 + Q] [(Q - l)g. + (3 - Q)g±(a)] g d nj D C 



\J [cos 2 (a) — sin 2 (a) cos 2 (0)] J 



^1 =V 2 [90 + 9s - g±(a)} n 3D C J — 
lj?=2\1 - Q 2 \\g^(a)\n 3D C 



nk 



(12) 

(13) 
(14) 



where a (8 )is the angle between B and n, (k and x). In 
the above, we write out the explicit a dependence of Eq. 
Has 

2-7T 

SJ-(a) = 9s - -j9d [1 - 3 cos 2 (a)] (15) 

The spin mode scales with y/k and is highly anisotropic 
due to the long-ranged and anisotropic nature of dipolar 
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FIG. 3: Collective mode spectrum for a 
and q — q c + Sq with q c = —0.87 and 8q - 
Hz from top to bottom. 



= 0.24tt, 6 = it/2 
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interactions. The spin mode can develop a roton mini- 
mum indicating a tendency towards crystalline order and 
has a strong dependence on a. The charge mode scales 
with k and describes phonon excitations of the supcrfluid. 
Notice the superfluid velocity depends primarily on the 
spin-independent contact interaction go. The magneti- 
zation mode is gapped and describes fluctuations in the 
magntiude of F. Notice the gap can vanish as a function 
of Q and a. 

We plot a representative collective mode spectrum 
in Fig. [3] illustrating the spin and charge modes for 
a = 0.247T, 6 = ir/2 and q near q c = -0.87 Hz. Notice 
the appearance of a roton minimum and the softening 
of the roton gap as q approaches q c . When q is below 
g c , the spin mode becomes imaginary at finite wavevec- 
tor indicating a dynamical instability towards periodic 
crystalline order. 

We then analyze the imaginary part of the collective 
mode spectrum in momentum space as in the left of Fig. 
21 This indicates the structure of instabilities and we 
construct the phase diagrams for the spin and charge 
mode in Fig. 21 The uniform ferromagnet mean-field 
solution exists within the green lines and the polar state 
with no magnetization is outside. Dynamical instabilities 
above the red line imply the mean-field state is unstable 
to small fluctuations. In general, spin mode instabilities 
are stronger with larger Im than charge mode instabil- 
ities except near the boundaries to the polar state where 
they are comparable. 

Notice the boundaries to the polar state intersect near 
a c = 0.357T indicating the uniform ferromagnet is not 
stable for any q even at the mean-field level. This is anal- 
gous to the magic angle effect familiar from NMR [26j |. 
At this magic angle, g±{a c ) — and rapid precession of 
the spins cancels out the combined effect of dipolar and 
spin dependent contact interactions which stabilize the 
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FIG. 4: Spin mode in momentum space (left) illustrating regions of unstable modes with shading indicating magnitude of Im uj^, 

and B the magnetic field for current experimental parameters a = it/2, q = 23]. The red marker indicates these parameters 
in the spin and charge mode phase diagrams (middle and left) which have similar plots illustrating regions of unstable modes 
in various phases (see text for description). The blue arrow indicates parameters for Fig. [3] 



ferromagnet. 

In the phase diagrams, D a (C a ) denotes regions of un- 
stable modes disconnected from (connected to) the origin 
k = and R indicates a ring of unstable modes. There 
are several trends to notice in the phase digrams. The 
uniform ferromagnet is stable only for a < 7r/4 away from 
the negative q boundary. For a < 7r/4 near the negative 
q boundary the spin mode is unstable along the direction 
perpendicular to B. This suggests an instability towards 
a striped phase modulated along this direction. From 
Fig. HI we see it is driven by transverse fluctuations SF 
of the magnetization F which modulation perpendicular 
to B. There is also an area of ring instabilities near the 
lower left hand corner. Generally, modes perpendicular 
to B are more unstable than modes parallel to B. 

For 7r/4 < a < a c , the primary instability is in the spin 
mode although the charge mode is also unstable in lim- 
ited regions near the polar state. Again, transverse fluc- 
tuations SF drive the spin mode instability towards mod- 
ulation perpendicular to B near the negative q boundary 
to the polar state. However, longitudinal fluctuations SF 
drive the spin mode instability towards modulation par- 
allel to B near the positive q boundary. This suggests 
striped phases perpendicular and parallel to B near the 
negative and positive q boundaries, respectively. In be- 
tween these boundaries, both transverse and longitudinal 
fluctuations SF are important. This indicates a tendency 
towards a checkerboard phase modulated along both di- 
rections. Above the magic angle a c < a, both types of 
instabilities are important throughout indicating a ten- 
dency towards checkerboard phases. 

The above analysis for fg, \k\ = also holds qualita- 
tively for general parameters fg, \k\ ^ 0. As an illus- 



tration, we plot in Fig. [5] the imaginary part of for 
the experimentally relevant parameters a = it/2, q = 0, 
/3 = in four cases: fg = and \n\ — (top left), 
fg = 0.8 and = K max (bottom left), /b = and 
|«| = (top right), fg = 0.8 and |«| = K max (bottom 
right). Here n max = 27r/60 /im -1 is the largest spiral 
wavevector (tightest winding) obtainable in current ex- 
periments [231 ] - In general, increasing fg suppreses the 
dynamical instability by decreasing of Im u)?. Increas- 
ing | ac | introduces additional instabilities with a length 
scale set by 2vr/| ac| which is typically 60 ^m or greater. 
They are well-separated from instabilities due to dipolar 
interactions which are typically at a length scale of 10 
/im. Experimentally, the instability also appears to oc- 
cur on shorter or longer timescales depending on whether 
the length scale of the spiral winding is smaller or larger, 
repspectively, than the largest length scale of the con- 
densate. We do not find a strong dependence of the in- 
stability timescale on the wavelength of the spiral wind- 
ing which was reported in [23| . One possibility for this 
discrepancy is that in experiments, spiral states are pre- 
pared dynamically by applying a strong gradient of the 
magnetic field. Winding the spins dynamically may in- 
troduce noise and fluctuations into the system and facil- 
itate the development of fragementation. 



DISCUSSION 

Our analysis shows the intriguing complexity of excita- 
tion spectrum of spinor condensates when dipolar inter- 
actions are taken into account. In particular we demon- 
strated that softening of the roton gap is associated with 
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FIG. 5: Imaginary part of the spin mode in momentum space 
with and without spiral order (right, left) as well as with and 
without longitudinal magnetization (bottom, top). Here we 
take current experimental parameters of a = tt/2, q = Hz, 
(3 = 0@. 



transverse fluctuations of the magnetization This lowers 
the dipolar interaction energy and such fluctuations are 
distinct from phonon excitations of the superfluid. Ear- 
lier studies on dipolar effects in polar molecules trea ting 
the spin degrees of freedom as frozen 

EE EE EE EE El 

obtained an excitation spectrum with one branch mixing 
the phonon and roton parts. 

Dynamical spin degrees of freedom have been con- 
sidered before in studies of ground state properties of 
spinor condensates but only with the bare dipolar inter- 
actions. Our analysis demonstrates rapid Larmor preces- 
sion and reduced dimensionality significantly modify the 
effective dipolar interaction. Using this effective interac- 
tion, we have shown uniform mean-field condensates have 
instabilities towards a possible supersolid with stripe or 
checkerboard crystalline order. 

Our work opens several new directions in the study of 
spinor condensates. We demonstrated that the uniform 
phase is unstable toward spin modulation for a whole 
range of wavevectors. However competition between un- 
stable modes and the resulting stable state still need to be 
explored. Possible candidate phases include both states 
with small modulation and states with spiral spin wind- 
ing. Another interesting question to be explored is the 
interplay of Bose condensation and spin modulation in 
the presense of thermal fluctuations. At the mean-field 
level spin modulation should appear simultaneously with 



the appearance of the condensate. However spin mod- 
ulation breaks translational symmetry and is a distinct 
symmetry breaking from the Bose condensation. Hence 
apriori the two may have different transition tempera- 
tures. 

The main significance of our work is developing a mi- 
croscopic theory which successfully explains the checker- 
board pattern observed in recent experiments of the 
Berkeley group [23] . The Fourier transform of the magne- 
tization observed in experiments had a distinct cross-like 
structure with a lengthscale of 10 fim which is in excel- 
lent agreement with the results shown in Figs. [4] and [5] 
Previous theoretical studies 27 , 2|| have addressed these 
experiments without taking into account the precession 
averaged dipolar interaction and have not explained the 
observed pattern of spin fragmentation. Although we 
analyzed the simplest geometry of an infinite two di- 
mensional layer, our results are in good agreement with 
the current experimental data. We expect however that 
the two dimensional model may not capture some effects 
which may be present in real systems such as pinning of 
the instabilities by the in plane trapping potential. 

In conclusion, we have considered the combined effects 
of dynamical spin degrees of freedom, Larmor preces- 
sion, and reduced dimensionality on the collective mode 
spectrum of spinor condensates with dipolar interactions. 
Starting from a mean-field state with uniform or spi- 
ral magnetization, we demonstrated the presence of a 
phonon and roton spectrum of collective modes. Soften- 
ing of the roton gap suggest instabilities towards a possi- 
ble supersolid state with stripe or checkerboard periodic 
crystalline order. 



METHODS 



Dipolar Interaction 

Here we outline the derivation of the effective two- 
dimensional time-averaged dipolar interaction. The ef- 
fect of the transformation in Eq. 2] is clearest in the 
bosonic coherent state path integral formalism for hi the 
evolution operator 

W = y V^V^e^ dtiC , C = i J cfx&dtV-H (16) 

where TL is given in Eq. [TJ Only the time derivative, spa- 
tial derivative and dipolar interaction tensor transform 
non-trivially. The transformed time and spatial deriva- 
tives 



dt -» d t + B B ■ F, 

V 2 V 2 .\k\k- V 
•— — 



B ■ F 



B ■ F 



(17) 



2m to 2to 

can be absorbed into Hamiltonian as modified linear and 
quadratic Zeeman shifts. 
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For the dipolar interaction tensor, the transformation 
needs to be expressed as one acting on K l ^ D (x) instead 
of ^3. Recall Fj k are 5*0(3) Clebsch- Gordon coefficients 
projecting the tensor product of two spin-1 representa- 
tions (lower indices) onto the spin-1 component (upper 
index). In particular, 



(18) 



for an arbitrary 5*0(3) rotation R implying R acting on 
the lower two indices is equivalent to R acting on the one 
upper index. This gives 



h^ D (x x ) — 



-tt/Bo 



Bndt 



tt/Bo 



R(t,S)f i ,h^(x-x')R(t,x') fj 
(19) 

where the bar denotes time-averaging and the explicit 
dependence on t and x + x! is removed as a result. Here 
\i^ D (x) is given in Eq. O The short-distance singularity 
I a?| — 3 in the Fourier transform 



and the following function 

gSdxdx>e-^-*\f{x)f{x') 
W{q) ~ 2jdxP(x) (28) 

the only quantity dependent on f(x). The asymptotic 
behavior of w(q) is 



"(9) = \ l^lj^l 



1 ~ £±l q 2 q<1 



D, 



-Q 



q > 1 



(29) 



with D = J dxf 2 (x), D +n = J dxdx' f(x)f{x')\x - x'\ n , 
D- n = J dxf(x)(—id x ) n f(x). For a gaussian trial wave- 
function 



fix) 



-x 2 /2 

/2n ' 



){q) = 2q / dke-( k2+2kq ^ (30) 
Jo 



and we stress the qualitative behavior of w(q) is rather 
insensitive to the detailed form of f(x). 



h l i D {k) = J d 3 xe-^ D (x) 
is regularized with the prescription 




(20) 



(21) 



and taking b — > at the end. 

To take into account confinement along n to thickness 
d n , we consider a general trial wavefunction of the form 
in Eq. [3] with 



P(*n) = -ff \ ^J- 



(22) 



where x n is the coordinate along h and f{x) is normalized 
to / dxf(x) = 1. This results in 



J dk n dxdx'h^ D (k, k n )e lk ^ x ~ x ^f{x)f{x') 
JdxP(x) 



(23) 

with k (k n ) the two-dimensional (one-dimensional) coor- 
dinate perpendicular (parallel) to h. We obtain Eq. [5] 
with 



h B D (k) = l-3g(k,0) 



(24) 



^ + - A g{k,+\K\k) + -^g{k,-\n\k) (25) 



_, 3 3 
h 2D( k ) = -^9(k,+\K\k) + -g{k,-\n\k) 

and the function g(u, v) given by 

2 i^\u+v\d n 



(26) 



g(u, v) 



D 



U + V 



w(\u + v\d n ) 



u + v\ 

! ( B ■ h) [1 - w(\u + v\d n )) 



dq 



Collective Modes 

Next we outline the mean-field and collective mode 
analysis. Taking — *S? with all quantities independent 
of x in Eq. [7] and extremize the free energy at fixed 
density n and longitudinal magnetization fg. The free 
energy is independent of 77 and 4> and we find two classes 
of solutions 



u = 0, 

IT 

2' 



v 



)t 3 + (1 - Q)r = f E 
+ (1 - Qlr- 1 = f E 



(31) 
(32) 



with sin(p) 2 = (r + r 



The constraints • 



sin(p) 2 < 1 and — 1 < r < 1 select the unique root for 
the cubic equations 



2Qsinh 

At 

IB 



rarcsinh 



fl /3 

2Qsign(/ B ) cosh 
2Qsign(/ s ) cosh 



(20^) 



Iarccosh(^ 



3 arCC ° Sh V2QQ3 



(\h 



Q < l 

= 1 
> 1 
(33) 

(34) 



w(q) 



\u\d„ 



(27) 



with Q = - 0I/3Q and Q given by Eq. E for the 
v = 0, 7r/2 solutions, respectively. For g±(a) > or 
a > a c (g±(a) < or a < a c ), the v = (y = n/2) 
solution has a lower mean-field free energy. However, we 
stress both solutions are stationary points for all a. 
The above constraints also imply the v = solution 

exists for Q > —\ (l + y/1 — f'A while the v = it/2 

solution exists for Q < — | ^1 — \J\ — f 2 ^j . In general, 
both types of solutions support transverse magnetization 



8 



except the fully polarized state along fg = ±1 and the 
polar state along Q > 1, /g = (Q < 0, fg = 0) for 
^ = (y = 7r/2). We focus on the r\ = solution for —1 < 
Q < +1 as it support a finite transverse magnetization 
at /s — as observed in experiment. 

Turning to the collective mode analysis, the linearized 
equations of motion are given by Eq. [9] with 



K 2m 



• k 



BF 



Q ' 



BF 



90 



in i \ 2m , 

p(0)C&* + g o p(0)C^ + 2.g s p(0)C*** T 



g d p(0)Ch% D (0)Tr 



TV- = 5oP (0)C** T + «? s p(0)C* T * 



(35) 



(36) 



with C = l/y/2 and h£ D (k) by Eq. 
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